Mutual friction in a cold color flavor locked superfluid and r-mode instabilities in 

compact stars 
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Dissipative processes acting in rotating neutron stars are essential in preventing the growth of 
the r-mode instability. We estimate the damping time of r-modes of an hypothetical compact quark 
star made up by color flavor locked quark matter at a temperature T < 0.01 MeV. The dissipation 
that we consider is due to the the mutual friction force between the normal and the superfluid 
component arising from the elastic scattering of phonons with quantized vortices. This process is 
the dominant one for temperatures T < 0.01 MeV where the mean free path of phonons due to their 
self-interactions is larger than the radius of the star. We find that r-modes oscillations are efficiently 
damped by this mechanism for pulsars rotating at frequencies of the order of 1 Hz at most. Our 
analysis rules out the possibility that cold pulsars rotating at higher frequencies are entirely made 
up by color flavor locked quark matter. 
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Introduction. — According to Quantum Chromody- 
namics (QCD), at asymptotically high densities strong 
interacting matter cannot be described as made up of nu- 
cleons; at sufficiently high densities nucleons are crushed 
and their quark matter content should be deconfined [ij. 
It remains an open question whether deconfined quark 
matter can actually exist in neutron stars that have cen- 
tral density at most ~ 10 times that of ordinary nuclear 
matter or whether stars entirely consisting of quark mat- 
ter can exist. The latter case may be realized if strange 
quark matter is absolutely stable with a lower energy per 
baryon than nuclear matter. Such a quark star would 
consist of quark matter fluid all the way to the surface. 

If deconfined quark matter exists in compact stars it 
is likely to be in one of the possible color superconduct- 
ing phases 0. Neutron stars that are more than a few 
seconds old have temperatures of the order of tens of 
keV at most and such temperatures are smaller than the 
critical temperatures of color superconductors which are 
generically of order tens of MeV 0] ■ The color supercon- 
ducting phase that is energetically favored at extremely 
high densities and low temperatures is the color flavor 
locked phase (CFL) [s]. In this phase u, d and s quarks 
pair forming a quark condensate that is antisymmetric in 
color and flavor indices. The order parameter breaks the 
baryonic number U{1)b symmetry spontaneously, and 
therefore CFL quark matter is a superfluid. 

In this Letter we analyze the r-mode instability of com- 
pact stars entirely made up of CFL quark matter. R- 
modes are non-radial oscillations of the star with the 
Coriolis force acting as the restoring force, they pro- 
vide a severe limitation on the star's rotation frequency 
through coupling to gravitational radiation (GR) emis- 
sion '4'|. When dissipative phenomena damp these r- 



modes the star can rotate without losing angular mo- 
mentum to GR. If dissipative phenomena are not strong 
enough, these oscillations will grow exponentially and 
the star will keep slowing down until some dissipation 
mechanism can damp the r-modes. Therefore, the study 
of r-modes is useful in constraining the stellar structure 
and can be used to rule out some color superconducting 
phases Q. 

We consider only the low temperature regime (T ^ 
0.01 MeV) where the thermal and dynamical properties 
of the system are dominated by the contribution of the 
superfluid phonons ip associated with the spontaneous 
breaking of [/(1)b- In order to provide an efficient mech- 
anism for damping of r-modes in Ref. [G] the processes 
if <-> ifip and ifip <-> ifip were studied. The latter process 
is largely suppressed, while the former can be seen to pro- 
duce a mean free path that scales in terms of the quark 
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chemical potential, /i, and temperature as ^ '--^ For 
low temperatures the authors of 6] find that i exceeds 
the radius of the neutron star (of order 10 Km) and there- 
fore this process does not provide an efficient dissipative 
mechanism and cannot be responsible for the damping 
of r-modes. Actually in this temperature regime one can 
treat the phonon system as an ideal bosonic gas. 

In the present paper we study a different dissipative 
mechanism which operates in rotating superfluids due to 
the presence of quantized vortices. In this case a mutual 
friction between the superfluid and the normal compo- 
nent arises because of the scattering of phonons off super- 
fluid vortices. We estimate the time scale for the dissipa- 
tion of energy associated with this process and conclude 
that for low temperatures and for reasonable values of 
the vortex mass this mechanism is efficient in damping r- 
mode instabilities for stars rotating at frequencies v 1 
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Hz. On the other hand, cold pulsars rotating at higher 
frequencies are unlikely to be CFL quark stars. As a mat- 
ter of fact in the temperature range considered here the 
populations of electrons f?\ and gapped quarks 0, Q are 
exponentially suppressed and they cannot significantly 
contribute to dissipative processes. 

We leave for a future project the study of the effect of 
mutual friction in the dissipation of r-modes in a higher 
T regime of the CFL phase. An incomplete analysis of 
the damping mechanisms as arising from other dissipative 
channels has been carried out in Refs. 0, 

Throughout, we use natural units h, = c = Kb = 1 
and Minkowski metric 77^^ = (1,-1,— 1,-1). 

Vortex dynamics. — We review the description of the 
forces acting on a vortex line pioneered by Hall and Vi- 
nen ^E^. 

The relative motion of the superfluid component with 
respect to the vortex produces a Magnus force per unit 
vortex length given by 



F*^ = upsi-Vs - Vl) X z. 



(1) 



where k is the circulation, z indicates the direction of 
the vortex line, ps is the superfluid density, Vc, is the 
superfluid (local) velocity and is the vortex velocity. 

The normal component reacts to the vortex motion 
through a force which can be decomposed into compo- 
nents parallel and perpendicular to the relative motion 



F^ = D(v„ - vl) + D'z X (v„ - vl) 



(2) 



where v„ is the velocity of the normal component and 
D and D' are two coefficients that depend on the micro- 
scopic physics of the system. At low temperature, these 
coefficients can be computed by analyzing the scattering 
of quasiparticles off a vortex [lO| . 

At equilibrium these two forces add up to zero 
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0, 



(3) 



and from this equation one can determine the vortex ve- 
locity to be given by 



Vl = + a'(v„ - Vs) + az X (v„ - v^ 



where 



d2+(l-d^)2 ' 



1 — a =^ 



:fi + (i-di)2 



(4) 



(5) 



where d|| = D/ nps , and d^ = D' /nps- 

Since F*^ is the force acting on the vortex generated 
by the superfluid component, the corresponding force on 
the superfluid is given by — F*^ and we will assume that 
there are no other forces acting on the superflow. 

Upon substituting Eq. (U) in Eq. ([T]) one finds that the 
mutual friction force on a unit volume of the superfluid 
component per unit vortex length is given by, 



^SN 



-iV„KPsQ!(vs-v„)-|-7Vi,Kpsa'[z;x (vs-v„)] , (6) 



where we have multiplied the force by the surface density 
of vortices Ny. This force is acting on the superfluid com- 
ponent due to the interaction of the normal fluid compo- 
nent with the vortices. The first term on the right hand 
side is responsible for dissipation. 

The expression in Eq. ^ is valid if the Magnus force 
and the mutual friction force equilibrate, i.e. if Eq. Q is 
satisfied. If we introduce a perturbation of the superfluid 
velocity it is not guaranteed that these two forces will 
equilibrate and we will have that the force per unit vortex 
length acting on the vortex is given by 

d5^fL 



SY" EE ,5F^ 



dt 



(7) 



where (Svl is the variation of the vortex velocity induced 
by the perturbation of the superfluid velocity and is 
the mass per unit length of the vortex (we will loosely 
refer to this as the vortex mass) . 

We will assume that when an external perturbation is 
applied to the system the normal component remains in 
thermal equilibrium. On the other hand, the superfluid 
and vortex velocities will have an oscillatory behavior 



(8) 



where lu is the frequency of the r-mode perturbation, 
which is of the order of the rotation frequency of the 
star 0]. Henceforth we shall drop out the "(0)" from 
the subscripts. We find that the force acting on a unit 
volume of the superfiuid component becomes, 

g-pSN ^ NyKpsi-aSvs + a'z x (5v,) , (9) 
where a', a are defined as a', a with the replacement 
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Phonon-vortex scattering. — The rotating CFL super- 
fluid is threaded with vortex lines whose density per unit 
area is Ny = ^ , where f2 is the rotation frequency of the 
star and where the quantized circulation is n 
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In order to obtain the coefficients D and D' appearing 
in the friction force in Eq. ^ for the cold CFL matter 
one has to evaluate the cross section associated with the 
scattering of phonons with vortices. We will restrict our 
analysis to the elastic scattering because inelastic pro- 
cesses are suppressed, see and consider phonons as a 
bosonic ideal gas. In this case one can describe the prob- 
lem of phonon-vortex interactions using gravity analogue 
models [H, Q EE | . 

The action of a phonon field moving on the top of a su- 
perfiuid background is the same as that of a boson prop- 
agating in a given curved space-time. One can consider 
the superfiuid phonon as propagating along the geodesies 
of the metric 



ds^ 



Qaijdx^dx'^ 



where Q^^ is the relativistic acoustic metric 15 1 



0/^1/ = Tlf^f + - 1) v^v^ , 



(10) 



(11) 
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which depends on the speed of sound Cs and the super- 
fluid velocity Vfj,. In cylindrical coordinates {t, r, 9, z) far 
away from the vortex one has that ~ (1,0, 0) and 
one can derive the cross section for the phonon-vortex 
scattering from the equation of motion of the phonon field 
in the corresponding curved background. For a phonon 
of energy E the differential cross section per unit vortex 
length for the phonon-vortex scattering is given by 



da Cs cos 6 , 9 ttE 



de 2nE tan2 f A 



(12) 



where 1/A = (l - c^) K/(27rc^). Notice that for the CFL 
superfluid there is an extra factor (1 — c^) in the expres- 
sion of 1/A with respect to the result of 13, Uj due to 
the fact that the speed of sound is relativistic. 

The parallel and transverse cross sections for the 
phonon-vortex interaction, for E' <C A, are given by 



1 



, 2 K 



8c2 



= l-ct 



(13) 



where k is the phonon momentum. A similar re- 
sult has been derived for non-relativistic superfluids by 
Volovik [l3 | employing the analogy with cosmic strings in 
(3-|-l)d and also with the problem of a spinning particle 
in (2+l)d gravity. However, in order to derive the exact 
phonon-vortex cross section one has to consider correc- 
tions due to the short distance behavior of the velocity 
field JJ]. 

From the expressions of the cross sections in Eqs. 
one can determine the coefficients D and D' that appear 
in the mutual friction force (see e.g. [IG]) 



£) = cri 



^sPn 



D' = aj_CsPn ; 



(14) 



where p„ is the phonon density. 

Mutual friction time scale. - We now discuss the dissi- 
pation of r-modes due to the mutual friction force in the 
CFL phase at small temperatures. In the reference frame 
where the normal component velocity is zero the friction 
force due to a superfluid velocity perturbation is 



(15) 



and the corresponding energy density dissipation is 



d£ 



5Fi^, • = ~N,,kpM5ws ■ Svl). (16) 
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The typical dissipation time is given by 

2nRe{a) , (17) 



1 _ 1 /d£ 

TMF 2E V dt 
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where we used the fact that £ = ^ps{6ws ■ Sw*). Since 
the expression of a depends on the ratio x = we will 
consider different values of this quantity. 



In case of negligible values of the vortex mass, i.e. x 
1, we have that 



27r5 



a ' 



405c6 



T 



(18) 



and taking Cs — 1/ \/3 wc obtain 



^— = 2ar2 ~ 2 X 18.1 f-^ n. (19) 

TMF V M / 

Since T <^ p, the dissipation time scale in this case is 
much longer than the period of rotation of the star. If 
we take the opposite limit, a; 3> 1, we have that ~ 
afl/x'^ and the corresponding time scale is even longer. 
The only case where dissipation time scale is of the same 
order or less than the period of rotation of the star is for 
a; ~ 1. Estimates of the value of the vortex mass give 
niy ^ p^"^, where ^ is the radius of the vortex core. The 
value of ^ is uncertain [16], but considering the largest 
possible estimate where ^ ~ V^, with A > 10 MeV the 
superconducting gap, we obtain that x < 10^^^ and then 
the effect of the vortex mass should be negligible. 

Critical frequency — Now that we have evaluated the 
dissipative timescale for the mutual friction, we can cal- 
culate the critical frequency above which the r-mode in- 
stability appears. 

We need to compare the timescale for the mutual fric- 
tion with the "growth" timescale for the r-modes tcr- 
For a quark star with a polytropic equation of state with 
index n = 1 one has that [4] 



1 



1 



TGR 3.26 \t:GnP 



(20) 



The condition of stability is > which for 

Tmf tor ' 

X <^ 1 translates into 



1/ < zy, ~ 2 X lOM - I Hz . 



(21) 



for a canonical star with M — 1.4Mq and i? = 12 
Km. Gn — 6.6742 x 10~^ cm-^g^^s^^ is the gravitational 
constant, and v^. is the critical rotation frequency. For 
T = 0.01 MeV, p = 400 MeV, we get - 1 Hz. In the 
unlikely case that the vortex mass is such that a; ~ 1, we 
obtain Vc — 10® Hz. 

Discussion. — Mutual friction in the cold CFL super- 
fluid can damp the r-modes if the rotation frequency 
of the star is less or of the order of 1 Hz. R-modes in 
cold CFL stars rotating at larger frequencies cannot be 
damped by this mechanism. 

We have obtained this result studying the interaction 
arising from the elastic scattering of phonons off super- 
fluid vortices. In the regime T ^ 0.01 MeV mutual fric- 
tion forces due to other quasiparticles present in CFL 
quark matter are exponentially suppressed and do not 
contribute signiflcantly in damping r-modes oscillations. 
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Mutual friction in conventional neutron stars has been 
considered in Refs. 
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18l |. The analysis in that case is 
more complicated due to the existence of both neutron 
and proton superfluids and an entrainment effect between 
them. The relevant microscopic process is the scatter- 
ing of electrons off the magnetic fields entrapped in the 
core of the superfluid neutron vortices. In that case, the 
associated damping times depend on the poorly known 
entrainment parameter. In some range of parameters the 
resulting forces are strong enough to damp effectively the 
r-modes. 

Let us remark that a different possibility not ruled out 
by the present analysis is that cold pulsars rotating at 
frequencies larger than are stars that have only a CFL 
core and the r-mode instability is damped by processes 
arising in matter constituting the remaining part of the 
star. 

The majority of observed pulsars have considerably 
larger frequency than 1 Hz [l9|, which indicates that it 
is unlikely that these pulsars are entirely made of CFL 
quark matter. Less than 25% of the observed pulsars 
have rotation frequencies less than 1 Hz [l^l; being iso- 
lated pulsars, their masses are unknown. Other observ- 
ables should be studied in order to see whether these 
slowly rotating stars are potential candidates for con- 
taining CFL matter. 

An alternative picture is that cold rapidly rotating pul- 
sars are quark stars which do not consist of CFL quark 
matter. When color neutrality and the effect of the 
strange quark mass are taken into account a promis- 
ing candidate is the three-flavor crystalline color super- 
conducting (CSC) phase [23|. This color superconduct- 
ing phase has gapless quarks as well as electrons that 
can provide r-mode damping. Moreover CSC quark mat- 
ter, differently from CFL quark matter, has a large shear 
modulus [m which is a necessary ingredient in explaining 
pulsar glitches. Another possibility is that the CFL-iiT" 
phase is realized [2^] , however the contribution of kaons 
to dissipative processes becomes relevant for tempera- 
tures T ^ 1 MeV [2^. It would be interesting to study 
mutual friction in the cold regime of this phase because 
of the presence of other light degrees of freedom. 
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